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Abstract 



This paper constructs a kinematic basis for spin networks with pla- 
nar or cylindrical symmetry, by exploiting the fact that the basis 
elements are representations of an 0(3) subgroup of 0(4). The ac- 
tion of the volume operator on this basis gives a difference equation 
for the eigenvalues and eigenvectors of the volume operator. For 
basis elements of low spin, the difference equation can be solved 
readily on a computer, yielding the eigenvalues and eigenvectors. 
For higher spins, I solve for the eigenvalues using a WKBJ method. 
This paper considers only the case where the gravitational wave can 
have both polarizations. The single polarization case is considered 
in a separate paper. 

PACS categories: 04.60.Pp, 04.30.-w 



I Introduction 

This paper calculates the eigenvalues of the volume operator for 
spaces with planar or cylindrical symmetry. The volume operator 
plays a key role in Thiemann's construction of a finite Hamiltonian 
for spin network canonical quantum gravity PP 

The quantization of systems with planar or cylindrical symmetry 
has been studied using both local field theory [3 El HI E] and spin 
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networks [HI E]. Before one can construct a spin network Hamil- 
tonian for these systems, one must construct a volume operator. 
Bojowald has given a discussion of some of the difficulties involved 
in constructing this operator |S|. 

Systems with planar or cylindrical symmetry are the simplest 
where gravitational wave propagation can occur. Simpler systems, 
such as homogeneous cosmologies and black holes with spherical 
symmetry, have higher symmetry and fewer degrees of freedom, but 
do not allow gravitational waves. In a homogeneous cosmology jH], 
for example, every point is equivalent to every other point. One 
may choose any one point, or vertex, of the spin network as repre- 
sentative. The Hamiltonian changes intertwiners and holonomies at 
that one vertex only. Since the action of the Hamiltonian is purely 
local, there can be no propagation of gravitational waves from point 
to point. For black holes with spherical symmetry, with or without 
matter jHl UHl HH 1121 UHl HI] , spherical symmetry rules out gravity 
waves. 

The construction of a loop formalism typically proceeds in two 
steps. The states of the Hilbert space must satisfy seven constraints. 
In the first step one constructs a basis for the Hilbert space which 
satisfies the six constraints which are easiest to treat (Gauss and 
spatial diffeomorphism constraints). The hardest constraint, the 
scalar or Hamiltonian constraint, is left for the second step. The 
basis constructed in the first step is the kinematical basis, and the 
dot product constructed in this step is the kinematical dot product. 
This structure is enough to determine the eigenvalues of the volume 
operator. 

This paper completes the first step. It constructs the kinemati- 
cal basis and dot product, and writes down equations which allow 
the eigenvalues and eigenvectors of the volume operator to be con- 
structed numerically. The paper does not go on to the second step 
and construct a Hamiltonian. 

The rest of this paper proceeds as follows. Section two describes 
the topology of the spin network. Section three defines the volume 
operator and sets up the kinematical basis. Section four writes out 
the eigenvalue equation for the volume operator in this basis and 
derives an analytical solution for eigenvectors with zero eigenvalue. 
For non-zero eigenvalues, and for basis functions of low spin, the 
equation can be programmed readily and solved on a computer. 
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Section five derives a WKB J method to find approximate eigenvalues 
of the volume operator. This approach was intended for high spins, 
but WKB J surprises (as it often does). It yields accurate eigenvalues 
even for basis functions of relatively low spin. 

This paper does not impose the constraints which limit the grav- 
itational waves to a single polarization. In classical gravitational 
theory, and in quantum field theories of gravitation, the single po- 
larization case is easier. However, in the spin network case I found 
the single polarization case to be harder. I will leave a discussion of 
that case for a separate paper. 

II Topology of the Spin Network 

A system with planar or cylindrical symmetry has two commuting, 
spacelike Killing vectors. In the planar case, if one shifts to the 
coordinates x,y suggested by the Killing vectors , the seven con- 
straints simplify considerably, because derivatives with respect to x 
and y may be dropped. (For the cylindrical case, replace the pla- 
nar coordinates x,y,z by coordinates 0,z,r. For simplicity in what 
follows, I shall discuss primarily planar symmetry, and devote only 
an occasional remark to the cylindrical case.) In particular, the x,y 
diffeomorphism constraints, and the X,Y Gauss constraints simplify 
enough that they can be solved and eliminated from the theory. Fol- 
lowing Husain and Smolin, I fix these four constraints by imposing 
the four gauge fixing conditions [7] 

E| = E| = 0; 

= Wy=Q. (1) 

Lower case letters x,y,z,. . . denote global coordinates; upper case 
X,Y,Z,. . . denote local coordinates rotated by the Gauss constraints. 
Setting the four constraints equal to zero and solving yields four 
more equations, 

= Pi.y = 0; 
Af = Ar = 0. (2) 

Both the triad and connection arrays are now block diagonal. A 2x2 
block contains fields with indices x,y and X,Y; a 1x1 block contains 
the field with indices z,Z. (For cylindrical symmetry substitute </',z,r 
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and $,Z, R for x,y,z and X,Y,Z.) The SU(2) local gauge symmetry 
has been reduced to local U(l), although the theory still contains all 
three generators Si of su(2). In order to express the consequences 
of the U(l) symmetry as clearly as possible, I choose the usual ma- 
trix representation of su(2) where is diagonal. I suppose that 
the reduction in components, eqs. ^ and (j2I) has been carried out 
classically, and I now set up a spin network formalism to quantize 
the reduced theory. 

In classical general relativity, the theory allows one to choose 
virtually any coordinate system; but it is not always clear which 
coordinates lead to the simplest equations. The choice of coordinates 
is not obvious unless the system has a lot of symmetry. 

Similarly in the spin network approach, the choice for the topol- 
ogy for the network is not obvious unless the system has a lot of 
symmetry. I will arrive at the topology for the planar case via the 
group theoretical approach favored by practitioners in the field of 
quantum cosmology. This approach has been discussed extensively 
in the literature ^3 US] , and for the most part my discussion will 
be a summary of results. 

As a first step in constructing the symmetry reduced theory, 
one uses group theory to construct a connection which embodies 
the symmetry. Modern treatments manage to avoid anything so 
declasse as solving the Killing equations, but in the end results are 
the same. The group theory determines the number of independent 
components and their form. It also determines the support of the 
connection fields. (In effect this step has been carried out already at 
eqs. (H} and (0), where we determined the number of independent 
components, and determined their support to be the z axis.) 

However, group theory by itself does not supply the topology of 
the spin network. For example, in a homogeneous cosmology, every 
point of the space is equivalent to every other point. Not supris- 
ingly then, the group theory predicts support for the homogeneous 
connection is limited to a single point. We know that, in the full 
theory, each vertex is connected to some number of edges. In effect, 
the group theory supplies a vertex (the point) but no edges. Mo- 
tivated by the full theory, one introduces edges, and promotes the 
connections to holonomies integrated over these edges. 

In the cosmological case, one can justify the introduction of 
edges d posteriori by constructing a Bohr compactified Hilbert space 
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which has some very unusual and useful properties (T7j. But it is 
unlikely anyone would have thought of doing this, had they not 
been motivated by the existence of edges and holonomies in the full, 
non-symmetric, theory. 

Similarly, in the planar case, the group theory predicts that sup- 
port for the connections is the z axis, but produces no edges. Moti- 
vated by the full theory, one introduces vertices along this 
well as two edges radiating from each vertex in x and y directions. 
The connections are then promoted to holonomies integrated over 
the edges. 

Holonomies on the z axis look like holonomies in the full the- 
ory. Each holonomy is integrated from one vertex to the next. 
Holonomies on the x and y edges are treated differently. These 
edges are given the topology of a circle; the two ends of each edge 
are identified. 

This identification is as first sight somewhat puzzling. (Most of 
the rest of this section records my initial doubts, and eventual ac- 
ceptance of circular topology for the xy edges. Readers who are 
comfortable with this choice of topology may wish to skip to the 
last three paragraphs of this section.) One would like to think of 
the symmetric theory as a reduction of a full theory. That is, when- 
ever one has a symmetry, one expects to start from the action of 
the full Hamiltonian on the full space, and end with the action of 
a simpler Hamiltonian on a smaller space. The ultimate example 
of such a reduction is a homogeneous cosmology, where one starts 
from a Hamiltonian acting on every point, and ends with a simpler 
Hamiltonian acting on a space consisting of a single, representative 
point. 

It is therefore natural to visualize the full spin network as con- 
structed from a basic unit, or module, which is repeated over and 
over to create the full space. For the homogeneous cosmology the 
basic module would be a unit cell containing only a single vertex, 
plus three edges. For the edges one could choose (say) three edges 
extending in positive coordinate directions, with the midpoint of 
each edge at the vertex. The basic cell would then be six half-edges 
radiating from a central point. Repeating this cell over and over 
generates the full space. (If the system possesses additional sym- 
metries such as isometry, one may need fewer than three edges. A 
change in the number of edges does not affect the present argument, 
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and I will ignore this possibility.) 

To return to my point: the modular picture does not immediately 
suggest an Si topology for each of these edges. Each basic cell 
contains an intertwiner at the vertex, plus six holonomies associated 
with the half-edges. When the full space is obtained by multiplying 
these basic units together, one index on each holonomy is already 
contracted with an index on an intertwiner; as for the remaining 
index, one might expect to contract it with the corresponding index 
on a neighboring cell. It is not obvious one should contract this 
index with another index in the same cell, which is what happens 
when one identifies ends and imposes an 5*1 topology on each edge. 

If one does not identify ends and contract indices, each cell would 
have six "danghng" SU(2) indices, indices not contracted with any 
other index. Those indices can be ignored when Gauss-rotating the 
cell, however. The uncontracted indices may be viewed as merely 
an artifact of splitting the full space up into identical cells. 

However, the "dangling index" picture does not hold up well when 
we consider reduction from the full to the symmetry reduced theory; 
circular topology seems essential. In more detail (and continuing 
with the homogeneous example for simplicity) , label each cell by its 
vertex v and label each half edge e by an index j^, which denotes 
the total spin of the holonomy on that edge (if one is using 2je + 
1 dimensional irreducible representations of SU(2) on each edge); 
or je denotes the number of spin 1/2 holonomies (if one is using a 
product of spin 1/2 matrices on each edge). By homogeneity, the 
two halves of a given edge must have the same je- Then the spin 
network wavefunctional for the full space is 



ip^ denotes the wavefunctional of a single cell. The full space is 
a product of such cells (the 11); and one must sum over possible 
assignments of edge spins (the E). The subscript m is shorthand 
for the six SU(2) indices on ip^. These indices are contracted with 
corresponding indices on neighboring cells. 

Since every point is equivalent to every other, the modules must 
identical, initially, and the dynamics must keep them identical. Dy- 
namics consists of identical transformations applied to every cell. 
In this spirit, I write the Hamiltonian constraint as a sum of con- 




(3) 
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straints, one for each unit cell. 

H = SvH 



The constraint = Hip implies 



= SjeHvV'v(je)m[nvYv^v'(jc)] 



m 



(4) 



I would like this to reduce to 



= EjeH^VvOe) 



i.e. to a single term in H acting on a single cell of ifj. However, this 
will not happen, in general, because the m indices on the square 
brackets in eq. (j3]) have ranges which depend on the je- The square 
brackets therefore do not factor out of the sum over je- 

I can, of course, eliminate the sum over je by confining myself to a 
single set of je', but this is too restrictive. I am left with the standard 
choice in the literature: identify opposite ends of each edge; trace 
over indices at opposite ends. Each cell then has no SU(2) indices 
to contract with neighboring cells. The square brackets in eq. (@)) 
are now identical, and they factor out of the sum. Equivalently, the 
traces change ^Z', eq. Q to 



where now there are no SU(2) indices on ipy. Each cell ipv{je) is an 
SU(2) scalar. 

I treat the planar case in the same way as the homogeneous case 
just discussed. The basic module is now the z line, plus two x and 
y edges at each vertex, with midpoints located at the vertex. The 
transverse edges are given the topology. For example, if I write 
out only the x edges and x indices at a given vertex. 



Here exi/2, for example, is the half-edge entering the vertex with 
intertwiner I. This structure may be rewritten in a manner which 
hides the traced indices. 



Those indices usually play no dynamical role anyway, since the vol- 
ume operator acts at vertices, i.e. its action affects only the mj, mj 
indices. 




(5) 



V'v ■ ■ ■ h[Ax, exi/2]]-[j/^jnjIjYiijmf h[Ax, Gxf/'2]mi,m' 



V'v = ■ ■ -Imi.mfhfA 

X) 6xf]mf,: 



,mi 
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Loosely, one can describe the rewritten holonomy as a "loop" 
holonomy which leaves from and returns to the same vertex. I will 
sometimes use this way of speaking, but note the holonomy is inte- 
grated over a line, not over a loop in the xy plane. The holonomy 
depends only on the connection A^. 

The final modular structure allows communication in the longi- 
tudinal direction, but no communication in transverse directions. 
Gravitational waves can propagate only along z. 

Ill Spin Network Basis States 

In the previous section I proposed a spin network topology consisting 
of an infinite line dotted with vertices; at each vertex two edges 
extend in the positive x and y directions with ends identified. I 
now associate a holonomy with each segment of the line and each 
transverse edge. For example, for the transverse x edge. 



The index B ranges over X and Y only. The expression on the 
right is not a trace. The h[Arc, e^f] matrix has two SU(2) indices. 
Both are contracted with corresponding indices on the intertwiner 
at the vertex; I am using the " loop" viewpoint described in the last 
three paragraphs of the previous section. The holonomies along z 
edges are similar to those in eq. (0), except that the su(2) valued 
connection is A^Sz and the integration is from one vertex to the 
next. 

The next step is to determine the selection rule which the holonomies 
must obey, because they preserve the residual gauge symmetry U(l) 
at each vertex. U(l) transformations are generated by the surviving 
Gauss constraint. 



A typical vertex, located at coordinate z, will have two z-edge holonomies, 
one beginning and one ending at z. The constraint eq. (|7j) generates 
the infinitesimal transformation Af Af — d^A which leads to the 
following finite transformation of the z holonomies. 







h[A^,e^i]h[A;„e,f] 
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= exp{i J rriziA^dz) exp{i J rrizfA^dz) 

exp[-iA{z){m^f -m^i)]H[A^]. (8) 

Here I have replaced Sz — or m^^; there is no point to retaining 
the Sz, since U(l) does not mix different eigenvalues of Sz- rrizf 
and ruzi are integers or half-integers, f and i label the outgoing and 
incoming z edges, respectively beginning and ending at the vertex. 

The same vertex will have two transverse edge holonomies begin- 
ning at the vertex. For the transverse connection A^., the infinitesi- 
mal transformation is 

A^Sc — ^ A^Sc + [A^ScjiASz], 

which implies the following finite transformation, when the connec- 
tion is promoted to a holonomy 

{exp[-iASz]h[A^, e^f] exp[+iASz])mf'mf 
= exp[-iA{mff - mf)]h[A^,er,f]mf'mf, (9) 

and a similar equation for h[Ay, eyf]nf'nf- 

Collecting together the x, y, and z contributions, I find that a 
U(l) transformation multiplies a vertex by the overall phase 

exp[-i(mzf - m^i + m^ - mf + np - nf)A(z)]. (10) 

U(l) invariance requires the relation 

2F := TUf — mji + nj — Ufi 

= nizf-nizi (11) 

Note that the intertwiner at each vertex no longer has to be a 
product of SU(2) 3 J symbols. The intertwiner can be any product 
of Kronecker deltas assigning specific values to the m's and n's, 
provided the selection rule eq. (fTTj) is obeyed at the vertex. 

It is straightforward to concoct a kinematic Hilbert space for the 
z holonomies, eq. (jH)). The functions exp[i m 9] constitute a complete 
set of functions periodic on [0,47r], for m integer or half-integer. 

However, I cannot use the SU(2) Haar measure to formulate a dot 
product for the transverse edge holonomies. To see this, I construct 
the simplest edge holonomy, that given by the spin 1/2 representa- 
tion of h. Denote the axis of rotation by n. Since h must lie in the 
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XY plane, it has the form 



h = (cos/3,sin/3,0). (12) 
for some angle 13. Then the spin 1/2 holonomy becomes 
/i(V2) ^ exp[m-cT^/2] 

= cos(^/2) +m-asin(^/2) 

= exp[-m^(/3-7r/2)/2]exp[i(Ty^/2] 

X exp[+ia,(/3-7r/2)/2]. (13) 

On the last two lines I have written the usual Euler angle decom- 
position for this rotation. It is clear that the rotation with axis 
confined to the XY plane depends on only two Euler angles, rather 
than the generic three. Let W denote the 2j-|-l dimensional matrix 
representation of SU(2). Then integration over only two angles 

j{\^'%,y{\i%^sme de d/3 

guarantees only a — h — a' — h' but not the stronger constraints 
a' = a, b' = b; the latter constraints must be satisfied in order for 
the 6 integration to yield j' = j. 

So far, I have been trying to follow a logical order: first construct 
an orthonormal basis; later, investigate the action of the volume 
operator on this basis. Logical order has not revealed a suitable dot 
product for the loop holonomies. 

It turns out it is better to reverse logical order. Start by investi- 
gating the action of the volume operator (on the spin 1/2 holonomy 
constructed above). This action suggests a basic structure which can 
be used to build higher order polynomials in the matrix elements 
of the spin 1/2 holonomy, while maintaining simple behavior under 
the action of the volume operator. Finally, construct the kinematic 
dot product suggested by the structure of these polynomials. 

I will need to define the action of the volume operator (more 
precisely the square of the volume operator). This operator is the 
product 

(V3)' = ezBcE|E|E^. (14) 

Each E operator must be integrated over an area, in order to make 
the volume invariant under spatial diffeomorphisms. In particular 
must be integrated over an area in the yz plane. For clarity I have 
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suppressed these areas in eq. p4|) . However, I cannot ignore them 
completely. Their precise extent is needed in order to determine the 
location on the yz area where E% can find an to grasp. 

That location is easy to find in the context of the full theory. 
The field has its support on an x edge, whereas the field is 
confined to the yz plane. The yz plane intersects the x edge only at 
the vertex. Therefore acts where the x edge holonomy meets a 
vertex. In the reduced theory, this means E^ acts at the two ends 
of the holonomy, eq. ©, where the holonomy intersects the vertex. 

I choose the yz area narrow enough in the y direction that the 
area includes only one vertex. When the full theory reduces to the 
symmetric one, neighbors in the y direction disappear anyway. 

In the z direction, I could choose an area which overlaps two or 
more vertices. Then E^ could act on one vertex at one z, while the 
other two triads in the volume operator act on a different vertex 
at a different z. I assume the yz area may extend in the z direc- 
tion halfway to the next vertices, but not all the way, so that the 
volume operator can grasp lines exiting from only a single vertex. 
Equivalently, I assume all three triads grasp edges exiting from the 
same vertex. This assumption has the advantage of simplicity. It is 
also reasonable, since propagation (which demands operators that 
change more than one vertex) is associated with components of the 
Riemann tensor, not with the volume. 

The E operators in a canonical quantization act like functional 
derivatives with respect to the corresponding A operators. When 
acting on an edge holonomy, the E bring down an su(2) generator at 
each end of the holonomy (never in the middle, because the volume 
operator acts only where the holonomy meets a vertex). For example 
the X holonomy is replaced by the following anticommutator. 



Subscripts A,B = X, Y only. I have omitted the triple delta func- 
tions, because they are always cancelled by the area and line inte- 
grals associated with E^ and . The overall factor 1/2 is a relic 
of the integrals over the delta's, 1/2 because the deltas occur at 
the endpoints of the edge integration, jk, is the product of Immirzi 
parameter times Sir G ; h = c = 1. 




{l/2)[SAh[A„e,f] + h[A„e,f]SA)]-fK. (15) 
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The anticommutator means that the volume operator generates 
an infinitesimal 0(4) transformation. To see this, specialize eq. fjl5|) 
to the spin 1/2 case: h[Ax, exf]^h[Ax, e^f]^^^'^^ (see eq. (fT^ ): and 
5'A^o"yi/2. Now consider the following two sets of SU(2) transfor- 
mations, 

R:={U-'h[A,,e,ff/'^U}; 

B:={Uh[A,,e,ff/'^U}. (16) 

The R's are just ordinary rotations; their infinitesimal form is a 
commutator. The infinitesimal form of the B's, however, is an anti- 
commutator. The sets R and B are special cases of the set of trans- 
formations {U'h[A^,e^fY^/^^U, U' ^ U}. The set is SU(2){g)SU(2). 

Modulo fine points about covering groups, SU(2)0SU(2) is 0(3){2)0(3) = 
0(4). To exhibit the 0(4) structure, introduce the suggestive nota- 
tion 



ih[Ax,exf 



(17) 



(1/2) + 1X4, Xi — 1X2 

xi + ix2 — + ixi 

Since the matrices U' and U are unimodular, the transformation 
U'h[Ax^exf\'^^^'^^U preserves the determinant, which is 

/ 2 I 2 I 2 I 2\ 

-(X1 + X2 + X3 + X4). 

The transformations R turn out to be ordinary rotations; they 
leave x^ invariant. A transformation B which has axis along direc- 
tion i rotates Xj into X4 (1 < z < 3). (The B stands for boost. The 
B's are of course rotations, not Lorentz boosts, but they become 
Lorentz boosts when their angular parameter is continued to a pure 
imaginary value.) 

If possible, I would like to avoid using 0(4) spherical harmonics as 
a kinematic basis. 0(4) is double trouble: the harmonics are prod- 
ucts of two rotation matrices, and the Clebsch-Gordan coefficients 
are products of two SU(2) Clebsch-Gordan coefficients. Fortunately, 
it is possible to use harmonics of an 0(3) subgroup of 0(4). 

To identify the 0(3) subgroup, I calculate the components Xj of 
h[Ax,exf\'^^'^\ eq. (fTTj) . The following traces give components xi 
through X4. 



Xi+ix2 = Tr[ih[Ax,exf]^'^^'^'^a^ 

= -sin(^/2)exp[i/3]; (18) 
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IX, = Tr[ih[A,,e,ff^h]/2 

= icos{9/2); (19) 

XI-IX2 = Tr[h[A„e,ff/^^a_]/V2 

= -sm{9/2)exp[-ip]] (20) 

xs = Tr[ih[A„e,ff/^^a,]/2 

= 0. (21) 

I have introduced the usual raising and lowering operators 

a± = ((Tx ± i(Jy)/V2. (22) 

In computing the above traces, it is convenient to write eq. (|T7jl as 
X ■ a + 1x4. 

From eq. (pT|l . the component X3 vanishes. This is not an ac- 
cident. There is no cxs in the expansion of eq. (fT!?jl because the 
rotation is confined to the XY plane. Further, is the only Sz 
valued connection left in the theory, after the gauge fixing, eq. (j21); 
and neither the Gauss constraint nor the volume operator have the 
power to change A^ into an Ar^^y. Hence the relevant operators main- 
tain X3 = 0. This suggests we do not need the full 0(4), but only 
the little group which leaves invariant the vector 

(Xi,X2,X3,X4) = (0,0,1,0). (23) 

The following theorems will make this idea more precise. First 
I recall some standard 0(4) theory in order to identify the genera- 
tors of boosts and rotations. For the SU(2)0 SU(2) transformation 
U'h[Axexf]U, let s' and s denote the generators of U' and U respec- 
tively, so that 

r / /] • / 

Pi, SjJ ~ l^ijkSk, 

[si, Sj] = ieijkSk; 
K Sj] = 0. (24) 
Then the generators of boosts and rotations are given by 



Proof: write s and s' as 

s = 
s' = 



bi = isi + s'^/2; 
n = {si-s[)/2. 



(s + s')/2 + (s-s')/2; 
{s + s')/2-{s-s')/2. 



(25) 
(26) 
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Demand that the general transformation U'h[Ax,exf]U reduce to 
the special transformations defined in eq. ()16|). For the boosts, the 
matrices U' and U must have the same generator, which means 
only (s + s')/2 can contribute; and similarly for the rotations, only 
(s — s') /2 can contribute. 

Next I need to identify the little group. Let (r'i,r2,r3) be gener- 
ators of su(2) and let V transform like a vector under the rj, that 
is 

[ri,rj] = ieijkrk; [ri, Vj] = ie^kYk. 

Then 

(ri,r2,r3) = (61,62,^3); 

{VuV2,Vs) = (-X2,+a;i,X4). (27) 

That is, the generators (61, 62, ''"3) of o(4)have the same Lie algebra as 
0(3); and the quantities (— X2, +xi, X4) rotate like a vector under the 
action of these generators. Note that X3 does not rotate; {bi,b2,r3) 
generate the little group. 

On the second line of eq. (j27p . the subscripts 1 and 2 have been 
interchanged, equivalent to a ninety degree rotation of the Vs. The 
rotation is necessary because bi (for example) is not quite ri. bi 
rotates xi into X4, whereas ri rotates V2 into V3. A rotation is needed 
to exchange V2 and Vi before the isomorphism will work. For later 
convenience I rewrite this isomorphism in terms of eigenstates of the 
ninety degree rotation: 

(V±,V3) = (±ix±,X4). (28) 

Proof of eq. (j27|) : it is straightforward to verify that the (61, 62? ''"3) 
obey the Lie algebra of o(3), by using the definitions eq. of the 
b's and r's, together with the commutation relations eq. ()24|) for 
the s and s'. To verify that the x's rotate like a vector, note that 
the U in U'h[Ax,exf]U is exp[(i(6i, 62, ''"3) ■ a], while the U' is the 
same, except for — «3. Therefore for the bi, the infinitesimal 

transformation is an anticommutator (with o"j/2, since h[Ax,exf] is 
spin 1/2); for the infinitesimal transformation is a commutator. 
bi (for example) generates the following infinitesimal transformation 
5[ih[Ax,e 



(1/2). 



[ih[A.,e.f](i/2),iai/2] 



1X\ ^4 

— X4 ixi 



(29) 
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or 

(5xi, (5x2, (5x3, 6x4) = (-X4, 0, 0, +xi). (30) 
Compare this to [ivi, Vj]^ = —eijkVk, or 

5V2 = -V3; (5V3 = +V2. 

After relabeling these equations as required by eq. (|77|). they become 
eq. ()30|) . One proceeds in similar manner to prove the rest of the 
isomorphism. 

I will refer to the vector on the right in eq. (j^Hj) or eq. as the 
basic vector. By inserting eq. (fTHj) through eq. (pljl into eq. (j^ . 
one can express the basic vector in terms of the angles (3 and 6. 

{V±,Vs) ^ (±^X±,X4) 

= (sin(^/2)exp[±i(/5-7r/2)],cos(^/2)) (31) 

Evidently the basic vector is a unit vector. Also, we can con- 
struct Condon-Shortley spherical harmonics Yim from the compo- 
nents of the basic vector pl^, but these spherical harmonics will 
have non-standard arguments: Yim{9 /2, (3 — t^/'2) rather than the 
usual Yim(6', /?). 

The next step is to build up more complex holonomies at the ver- 
tex: multiply together L matrix elements of e^^/]^/^ to form all 
possible homogeneous polynomials of order L in the matrix elements. 
These polynomials form a rank L reducible representation of the 
little group 0(3). Since we are multiplying together L identical vec- 
tors, we must break into irreducible representations by symmetrizing 
and taking traces; antisymmetrization gives zero. The irreducible 
representations are then just the Y^'m with L' = L, L — 2, L — 4,. . . . 

Therefore the set {Yirni^/^, P — 7r/2)} is complete, and I adopt 
it as a basis. Because these have argument 6/2 rather than 6, 
the usual dot product for the spherical harmonics 

/ dp / sin e d9 
Jo Jo 

will have to be modified slightly. 

< I'm' I Lm >= r"d/5 r sm{e/2) d{e/2) Y*,^,YLm. (32) 
Jo Jo 

I should perhaps emphasize that I have actually constructed two 
bases and two dot products, one for x holonomies with angles 9^, (3^, 
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and another for y holonomies with angles 0y,(3y. These angles are 
unrelated except in the one-polarization case, where classically, the 
two axes of rotation are at right angles, /^j; = /3y±7r/2. In that case, 
the kinematic Hilbert space and dot product must be reconsidered 
from scratch. 

Note that the basis is far easier to use than (say) the set 
{V\—(3 + 7i/2,9,P — 7r/2)om}, which would be a straightforward 
generalization of eq. (fT^ to a rotation matrix of higher spin. The 
Pom orthonormal, but they are not representations of the 0(3) 
little group, and the set {I^om} is not closed under a grasp by the 
volume operator. 

As for the lim, a grasp by E^, for example, multiphes the basic 
vector by the spin 1 representation of 5*+; this basic map induces a 
map of the higher harmonics into themselves, a map given by the 
spin L representations of S+. Symbolically, 

E^Yl^ = S^,Yw < L, m' I S+ I L, m >, (33) 

with dHl 

< L,m± 1 I Sa = S± I L,m >= ^(L ^ m)(L ± m + l)/2. (34) 

I have now completed the construction of the kinematic basis and 
dot product; and described the action of the volume operator triads 
on this basis. 

IV Eigenvalue Equation for the Volume Opera- 
tor 

To summarize the results of the previous section: I now have a basis 
set of holonomies; at each vertex there are two z-edge holonomies 

H[Af] = exp[iy" m^iAfdz] exp[i^m2fAfdz] (35) 

plus two "loop" holonomies 

Yl.„..(^x/2,/3x - vr/2)YL,„,,(^y/2,/?y - n/2). (36) 

At eq. (fTTj) I worked out the conseqences of the surviving U(l) 
gauge invariance, but for standard holonomies hm'm rather than the 
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new Yirn basis. It is easy to see from eq. (fT^ that a gauge rota- 
tion changes only f3 while leaving 6 alone; therefore to discover the 
gauge behavior of the Y's we need to study their /? dependence. 
The holonomies h^'m have /? dependence exp[i(m — m')/3] (see, for 
example, eq. (|13p). Except for a normalization, Yi^rn^ is the rota- 
tion matrix Vq^^, which has P dependence exp[zm^/?] Therefore at 
eq. the differences m - m' should be replaced by m^; and there 
is a similar replacement for the y indices. The U(l) selection rule, 
eq. (fTT|) simplifies to 

2F = iTLx + my 

= m^f - m^i. (37) 

The z holonomies of eq. (jH)) are eigenfunctions of the factor in 
the volume operator; and the remaining [E] operators in the volume 
map the loop holonomies eq. (jH^ into themselves as at eq. (jMj) 
Therefore the volume operator will not change mjinzf, Lx,ot Ly. 
It will not change F= {rrix + my)/2, F for "fixed", because of the 
selection rule eq. (jHTjl . The volume operator can change the quantity 

D:= (m,-my)/2, 

D for "difference". Therefore an eigenfunction of the volume oper- 
ator will be the product of the two z-edge holonomies eq. (jHI) times 
a sum 

I A; LxLyF >= SDYL,mxYLymyC(mx, my) (38) 

For simplicity I have suppressed the L and F dependence of the c's. 
The z dependence of the volume operator eq. (jl4p acts on each of 
these basis elements as follows. 

{VsyHlAz] I X;L^LyF> 

= (7^/2) (m^j - mzf)H[Az] 

X {V2f \ X; L.LyF > . (39) 

iffAzjis defined at eq. (jH)); the constants 7^/2 are as at eq. (fT3j) . 
(V2)^ is the determinant of the 2x2 subblock. 

To complete the action of V3, eq. (jH^ . I must determine how V2 
acts on the expansion eq. (j38j) . I shift from X,Y to the combinations 
X ± iY in V2, to simplify later matrix elements. 

El := (%±zE?.)/v^; 
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{V2f = ezABE^A^l 



= i{ElEl-W_El). (40) 

The i and the minus sign come from e^zp± = ±i. Note how plus 
indices are always contracted with minus indices. For example, 

Ej^ax + E^cTy = E^(T_ + W_a+ (41) 



E^ is essentially the functional derivative with respect to A^, there- 
fore when it acts on the basic holonomies (eq. ()13|1 . and eq. ()18j) thru 
eq. (j20|) ) . the functional derivative replaces each holonomy by 
its anticommutator with a+. (Again, plus indices always pair with 
minus.) The basic holonomies transform as a spin 1 representation 
under this anticommutation. When acted on by E^, therefore, the 
Y's transform as the 2L+1 dimensional representation of S±. 

{V2f\\-L,LyF> = i^^Kjlf 

X Sc(ma;,mj^)[yL^„^+i < L^m^ + 1 | 5+ | L^m^ > 

yrriy-l < Ly^y ~ 1 I ^- I ^y^y 
- (x ^ y)] 

= {-iKl2f\\\-L^LyF> (42) 

On the last line I have assumed the state is an eigenstate of (V2)^. 
Multiplying both sides by iL^m^lLymy, using the kinematic dot 
product introduced at eq. (j32j) . and the matrix elements given at 
eq. (jSH), I find the following eigenvalue equation for A. 

2Xc{m^,my) = ig^L^.m^, Ly.rriy) c{m^ - l,my + I) 
- ig{Ly,my,L^,m^) c{m^ + l,my - 1); 

g{L^,m^,Ly,my) = {L^ - + 1)(L^. + m^,) 

X ^JiLy + my + l){Ly-my). (43) 

In terms of A the eigenvalues A3 of the original operator (Vs)^ 
are, from eqs. and (j^ . 

A3 = {^K/2)'{m,i - m,f)A (44) 

Because F= {rrix + iTLy)/2 is held fixed, only the quantity 

D := (m. - my)/2 (45) 
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is incremented in eq. f)43p . The equation is an ordinary difference 
equation masquerading as a partial difference equation. To exhibit 
the ordinary difference character, I make the following replacements 

= F + D- 
my = F — D; 
c{m^±l,myTl) = c{D±l). (46) 

With these replacements, eq. becomes 

2Xc{D) = i^(L,-F-D + l){L, + F + D) 

X ^{Ly + F-D + l){Ly -F + D)c{D- 1) 

- i^{L, + F + D + l){L,-F-D) 

X ^{Ly + F-D){Ly-F + D + l) c{D + l) (47) 

I have been unable to find a compact, analytic solution to the above 
equation. The next three paragraphs, which describe my efforts to 
find such a solution, are perhaps of interest only to readers who are 
already familiar with the standard literature on ordinary difference 
equations [12]. Some readers may wish to skip these paragraphs on 
first reading. 

To make contact with the standard literature, which treats pri- 
marily equations with rational coefficients, I need to get rid of the 
square roots . This is easily done by a change of dependent variable. 

c{D) = ^J{L^ + m^)\/{L^-m^)\ 

X ^iLy + my)\/{Ly-my)\ d{D) 

= ^J{L^ + F + D)\/{L,-F-D)\ 

X ^{Ly + F- D)\/{Ly -F + D)\ d{D) (48) 

If this is inserted into eq. and square roots cancelled, the re- 
sulting equation for d is 

2\d{D) = i{Ly + F~D + l){Ly-F + D) d{D-l) 

- i{L^ + F + D + 1){L^- F ~ D) d{D + 1). (49) 

The coefficients are now rational, but they are quadratic functions of 
the independent variable D. Solutions when the coefficients are linear 
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are already available in the literature, but for quadratic coefficients 
one must construct a solution. 

One approach is to assume the solution is a series of factorials 
d(D) =S„a„/ (-D — n)!. (Factorials play the same role in the theory 
of difference equations that powers do in the theory of differential 
equations.) The problem of determining the d's is turned into the 
problem of determining the coefficients a„; for an equation with 
quadratic coefficients the new problem generally is as hard, or harder 
than the original problem. 

One may also try Laplace's method: write d(D) as the integral 
transform of a kernel, then show that solving the original difference 
equation is equivalent to solving a differential equation obeyed by 
the kernel. The differential equation corresponding to eq. (j49|) is of 
a type unknown to me. It has four regular singular points. If one 
attempts to solve the differential equation with a series solution, the 
recurrence relation for the coefficients in the series is as hard to solve 
as the original difference equation, eq. (jlUj) . 

Although analytic solutions are hard to find, numerical solutions 
are easy to implement. For low values of the L's, one may write 
eq. (jTTj) as a matrix equation, M-c = Ac, and solve for the eigenvalues 
of M. For larger values of the L's, I develop a WKBJ technique in 
the next section. For the rest of this section, I will discuss exact 
symmetries of eq. (jTTj) . 

Eq. (I47|) possesses the following symmetry, which is easy to prove. 
Write eq. (jTTj) in a matrix notation, M-c = Ac, and let {c{D; A)} be 
the components of a vector c which satisfies eq. ()47|) with eigenvalue 
A. Then the vector with components {{~l)^c{D; X)} also satisfies 
the equation, with eigenvalue —A. It follows that the eigenvalues 
occur in pairs (A, —A) , except for possibly the zero eigenvalues. 

This is a good place to mention that the operator needed by 
Thiemann P for his construction of the spin network formalism is 
actually the absolute value, | (Vs)^ | , rather than (Vs)^ itself; hence 
the relevant eigenvalues are | A3 |, and the theorem just proven 
states that eigenvalues of | (Vs)^ | are at least doubly degenerate 
(except possibly the zero eigenvalues). 

It is also possible to prove that there is at most one zero eigen- 
value, for given values of Lx,Ly, and F, and to construct the zero 
eigenf unction explicitly. When A = 0, eq. (jTTj) collapses from a sec- 
ond order recurrence relation to a ffist order relation. (The relation 
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connects every other value of D, relating c(D+l) to c(D-l); it is first 
order, with increment 2 rather than 1.) Square c{D + l)/c{D — 1) 
to get rid of the square roots; one then has a first order equation for 
the squares of the c's, with rational coefficients. This is a standard 
form, with solution known up to a normalization constant N • 



for D = max D, max D - 2, max D - 4, . . ., min D; and c(D) = 
otherwise. 

Eq. (j5(J|) must satisfy the boundary conditions that c(D) vanishes 
outside the limits max D >D > min D. To find the limits on D, I note 
that (although the symmetry is U(l)) I am using a basis of SU(2) 
spherical harmonics. Therefore D is constrained by the SU(2) limits 
—Li < rrii < +Li . These limits may be turned into limits on F± D 
by using the definitions F = (mj, + 'my)/2 ,D = (m^, — my)/2. It is 
then straightforward to derive the following limits on D. 

max(-Lx - F, -Ly + F) < D < min(+Lx - F, +Ly + F) (51) 

From these limits, the denominator f(D) in eq. ()5()|1 becomes infinite 
for D = max D + 2, max D + 4,. . . and for D = min D - 2, min D 
- 4,. . ., which enforces the boundary condition. Note the numerator 
is finite at those points. 

Since the recurrence relation connects only every other value of 
D, one might suppose there is another zero eigenvalue, with c(D) 
given by eq. ()50p for D = max D -1, max D - 3,. . ., min D +1; and 
D = otherwise. However, this solution does not obey the boundary 
conditions; e.g. c(max D + 1) is non-zero. 

Note also that the boundary condition requires the series D = 
max D, max D - 2, ... to terminate at min D, rather than min D + 
1. If the series terminates at min D + 1, then from eq. (j^Uj) c(min 
D - 1) will be non-zero, which violates the boundary condition. It 
follows that min D and max D must differ by an even integer, and 
the total number of allowed values of D, max D - min D -|- 1, must 
be an odd integer. This result is consistent with our earlier result 




(50) 
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that non-zero eigenvalues always occur in pairs (A, — A). For given 
values of {L^, Ly, F) and therefore given values of (max D, min D), 
there will be one zero eigenvalue if the number of allowed values of 
D is odd; otherwise there will be no zero eigenvalues. 

The eq. (j47|) links c's all having the same value of the parameter 
F. Put another way, there is one set of equations eq. (jTTj) for each 
value of F. A symmetry relates the equations for F to the equa- 
tions for -F, however, so that there is no need to solve both sets of 
equations. Let {c(D;F)} be a solution to eq. (jTTj) with eigenvalue 
A; then {c(-D;F)} is a solution to the equations eq. (jTTj) with F— >• 
-F and eigenvalue —A. (Temporarily I have restored the suppressed 
F dependence of the c's, for clarity.) The proof is straightforward, 
because changing (D,F) to (-D,-F) in the coefficients of eq. (|T7|) in- 
terchanges the two terms on the left, therefore changes the sign of 
the left hand side. (The c's on the left must be relabeled correctly; 
for example, c(D+l) becomes c(-D-l), not c(-D+l).) 



In this section I use results from a previous paper on WKBJ solu- 
tions to recurrence equations f^H", and I obtain a WKBJ solution 
for the recurrence relation eq. (j47|) . Since derivations were given in 
the earlier paper, for the most part I shall avoid derivations and 
motivate results using physical arguments. However, the earlier pa- 
per applied the theory to the 6J symbols, where there is no need to 
quantize an eigenvalue, hence no need to derive quantization condi- 
tions, formulas analogous to / pdx = {n + l/2)h in standard quan- 
tum mechanics. Later in this section I include some detail from the 
previous paper, enough to extend the theory slightly and derive a 
quantization condition. 

It is perhaps not surprising that the recurrence relation eq. (I47|) 
has a solution of WKBJ type. 



A recurrence relation may be turned into a second order difference 
equation, since the c's may be replaced by central differences. 



V WKBJ 



c(D) = Aexp(iS). 



(52) 




c(D + 1) -c(D); 

c{D + 1) - 2c{D) + c{D - 1); 
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c{D + 1) = c{D) + S^c; 
c{D-l) = -c{D) -6^c + 6^c. 



(53) 



Differences are very close to derivatives, tlierefore recurrence rela- 
tions are very close to differential equations. 

For a WKBJ solution to be possible, the recurrence relation does 
not have to contain the small parameter h. For example, WKBJ 
has been applied to the classical equations describing waves moving 
through an inhomogeneous medium, 

d^^/dx^ + k(x)2^ = 0. (54) 

A necessary condition for validity of WKBJ in this classical context 
is small derivatives, d"'k/dx"' =order kL'"", where L is a (large) 
length characterizing the rate of variation of the dielectric constant. 
In the present application, derivatives are replaced by differences; 
the quantities analogous to k are the square root functions in the 
recurrence relation eq. (@7j); and the large parameter (s) L are the 
quantities L^. ± F,Ly ± F. There is no , but we are definitely in 
the limit of large quantum numbers. 

I now state the requirements for a second order recurrence rela- 
tion 

g_(D, L) e(D - 1) + go(D, L) e(D) + g+(D, L) e(D + 1) = (55) 

to have a WKBJ solution [20]. The parameters L must be large and 
the coefficients g must satisfy the following conditions. 

{9+-g-)/g- < order (56) 
go/g± < order unity; (57) 
S^'g/g < order (58) 

Consider first eqs. and (jFTf) . If these conditions are not 
satisfied initially, often a change of independent variable will lead to 
g's which satisfy these conditions. The coefficients in eq. (jT7|) do not 
satisfy eq. (j^ . but I can remedy this by changing the dependent 
variable, 

c(D) = (i)°e(D). (59) 
This brings eq. (jTrf) to the standard form eq. with 

g-iD) = ^{L^-F-D + l){L^ + F + D) 
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X ^{Ly + F -D + l){Ly-F + D)- 
g+{D) = g-{D + l)- 

go{D) = -2A (60) 

Next consider the third condition, eq. (j^H)) . This is the analog of 
the " small derivatives of k" requirement and is the most important 
condition. It implies that the coefficients behave like polynomials 
under differencing, rather than like sinusoids, say. 

To test whether eq. is satisfied, it is convenient to approxi- 
mate differences by derivatives. If derivatives are falling off as 1/L, 
then so are the differences. 

6^g ^ (dg/dD)(AD = 1) + order(d^g/dD2). (61) 

When a square root in eq. is differentiated with respect to D, 
in effect the square root is divided by factors of order (L±F) ± D. 
Therefore the difference is down by a factor of order 1/(L±F) (not 
order 1/D; D does not have to be large). 

The g± will not obey eq. (j58p (will not behave like polynomials 
under differencing) whenever D is near zeros of the square roots, 
where the g± are non analytic and differences ~ derivatives can be 
badly behaved. 

Where are these zeros, and what is their physical interpretation? 
At eq. (I^Uj) I derived limits on D using the SU(2) constraints —Li < 
rrii < +Li. 1 recall these limits here: 

max(-Lx - F, -Ly + F) < D < min(+Lx - F, +Ly + F) (62) 

From a comparison of eqs. (j6Up and (j62|) the zeros of the g's occur 
at points where D is approaching the limits imposed by SU(2). I 
will refer to the limits on D given in eq. (j62p as the SU(2) limits on 
D, and the corresponding zeros of the g's as SU(2) zeros. I cannot 
expect a WKBJ solution to work near the SU(2) limits, and I will 
have to check that these values of D occur inside classically forbidden 
regions, where the solution is negligible anyway. I therefore need to 
locate the turning points, values of D where the solution shifts from 
classically allowed (sinusoidal) to classically forbidden (exponential), 
in order to check that the SU(2) limits are outside the turning points. 

To find the turning points, I apply formulas derived in [20] • In 
classically allowed regions, where the solution is sinusoidal, the am- 
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plitude and phase in eq. are given by 

= const. [((7+ + ^?_-2A)(-(7+ + 3(7_ + 2A)]-i/2. ^gg) 



5^S = arccos[(^_ + 2A)/2^_] (64) 
= 2arcsin^[((7_ + (7+-2A)/(4(7_)]. (65) 

When using these formulas for initial orientation, it is permissible 
to replace g+ + g- = 2g^, because from eq. (jUnjl g+ — g- = S^g-, 
which is assumed <^ g±. With this replacement, the above formulas 
simplify to 

^ const, [{g^ - A')]-^/^ ^gg) 
6^S ^ arccos[A/^_] (67) 
^ 2arcsin^(((7_-A)/(2(7„). (68) 



For the qualitative discussions in this paper I shall use these ap- 
proximate formulas. However, for numerical work one should use 
the more exact formulas. Neglecting the difference between g_ and 
g^ is equivalent to setting L(L + 1) = L^, which introduces signifi- 
cant errors even when the L's are as large as 10. 

In order to visualize the classically allowed region and the SU(2) 
zeros, it is helpful to make a rough plot of the function g^ versus 
D. From eq. (jEDI), g- is a quartic, therefore the sketch looks like the 
usual Mexican hat potential, g"^ has four zeros at the four SU(2) 
limits, eq. (jHH). Because of the min and max in eq. (jU^ . the SU(2) 
allowed region for D is the segment on the D axis lying between the 
two zeros closest to the center of the hat. Now draw a horizontal line 
at a height A^ above the D axis. From eq. ()66|) . the turning points 
are the two values of D (closest to the center of the hat) where this 
line cuts the hat. As A — > 0, these turning points approach the 
SU(2) limits and the WKBJ approximation breaks down. 

There is also a problem as the A^ line approaches the top of 
the hat. Recall the physical interpretation of the quantum number 
n in the usual quantum mechanical WKBJ formula / pdx = (n + 
l/2)h. n counts the number of half wavelengths which fit between 
the two turning points. In our case, as A^ grows large, the two 
turning points coalesce, and n becomes small; consequently WKBJ 
will be inaccurate. The "large" quantum number is not A^, but A^ - 
order(L^), since g^ is order (L^) near the top of the hat. It is perhaps 
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best to think of the Mexican hat as an upside-down potentiaL The 
"well" of the potential is at the top of the hat. 

From eq. ()68|) . I must take A > 0, in order for the sine to have 
a zero at the same time as the cosine eq. (|67p becomes ±1. I seem 
to have lost the negative A eigenvalues. I can recover them if I 
use {—i)^ instead of (i)^ in eq. (|^. This yields eq. (jSSI), except 
qq — i> —Qq. From eq. (|59|) this replaces A by —A everywhere. I have 
rederived a theorem from the previous section: the amplitudes for 
A and -A differ by (-1)^. 

I must now derive a solution near turning points, since the eigen- 
values A are quantized by the usual requirement that the WKBJ 
solution connects to exponentially decaying solutions at left and 
right turning points. Reference [201 derives the necessary connec- 
tion formulas. I change independent variable, 

e(D) = Z(D)/Vgr (69) 

and insert this form into eqs. ()55|) and ()59|) . I find that Z obeys 

6^Z + 2(g„ - A)Z/g„ = 0. (70) 

In deriving eq. (|7U|) I have assumed that I am near — = 0, 
but far from the SU(2) zeros of so that I can neglect higher 
differences of the y/gl in eq. . 

From the Mexican hat plot, there are two turning points D< and 
Dy, to the right and left of the center of the hat. Consider first the 
smaller turning point D = D^. I expand 

= + {D - D<)Z/k<; 
l/fc< = [2dig,)/dD]iD^)/X. (71) 

I have assumed that the zero of the second term in eq. ()70p is lin- 
ear rather than quadratic. This means 1/A;< cannot vanish. l/fc< is 
essentially the slope of the Mexican hat potential at the smaller turn- 
ing point, so is indeed non-zero, and also positive. At the smaller 
turning point, and in the classically allowed region, 

z< := D-D^>0; 
1/A;< > 0. (72) 

/c< is order L, assuming that the A and in the definition of /c<, 
eq. (j7ip . are coefficients in the same recurrence relation, therefore 
are of the same order in L. 
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Eq. ()7ip is a recurrence relation for the Bessel and Neumann 
functions. The Neumann solution can be discarded, because it has 
exponentially diverging behavior in the unphysical region, and I get 

e(^) = I c< I J-^<+2fc<(2/i;<)/v/^ 

^ I c< I cos[-2^4/A;</3 + 7r/4] 

X l/v/[A:<^< - \zj2f]Tig_. (73) 

The comes from eq. ()69p. On the second line I have used the 
Debye asymptotic limit, rather than the usual Hankel limit, because 
both the index and the argument of the Bessel function are large 
PT] . The form given is valid for 2/c< ^ z^> ?>^2k^, z large but not 
too large. I am free to choose the overall phase of the wavefunction, 
and I have chosen the constant c< to be positive. 

Now consider the larger turning point, D = D>, and again ex- 
pand 

= 6^Z+{D-D>)Z/ky, 
l/fc> = [2dig^)/dD]iD^)/\. (74) 

The slope on the right side of the Mexican hat, 1/A;>, is now negative, 
and so is z> = — in the classically allowed region. Again, the 
solution is a Bessel function 

e{D) = c>J+,^^2k>{-2ky)/y/gZ 

^ c> cos[-2^(-z>)V(-A:>)/3 + 7r/4] 

X l/^[fc>^> - (^>/2)2]7r^7- (75) 

The constant c> may have either sign, since I have already used 
up all the phase arbitrariness in the solution when I chose c< to be 
positive. 

eqs. (fTSj) and (f73|) must be matched up with the WKBJ solution 
eqs. (jU^ and The difference equation eq. has the following 
exact solution for S. 

D 

x=D^ 

= S{D^) + 5:2arcsin[^((7_ - X)/2g_] (76) 
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- S{D^)+Y^^[2{g_-\)/g_] 

= S{D^) + j24^<- (77) 

The third hne is arcsinx = x, vahd near turning points. The last 
hne uses the expansion of eq. (ffT|l . 

I would like to replace the sum in eq. ((TTj) by an integral, then 
integrate to obtain S (equivalently, replace 

^ dS/dD = dS/dz, 

then integrate). It is not immediately clear I may do this, because 
the square root in eq. ()77|) is not polynomial-like. However, consider 
the ratio 

r(z) = S^^^ov^/(2v^/3), (78) 

The numerator is the sum on the last line of eq. (|77p ; the denom- 
inator is the approximation to the sum obtained by replacing the 
sum by an integral. For D — = z = 3, r = 1.2; for z = 10, r = 
1.07. In words: as I move away from the turning point, the integral 
is dominated by regions where the function is polynomial-like, and 
the integral becomes a better approximation to the sum. Therefore 
replacing sum by integral in eq. (j77p is valid near z = D — (but 
not too near). I integrate eq. (f77|) to get S and then substitute into 
eq. (jSai: 

ReAexp{iS) = A cos[ T ^^Jk^ + S{D<)] 

= Acos[2y^z3/fc</3 + 5(D<)]. (79) 
Comparing this to eq. (f7H|l I get 

S(D<) = -7r/4. (80) 

Now consider the WKBJ solution near the larger turning point 
D = D>. 

D 

ReAexp{iS) = A cos[^(2 arcsin) — 7r/4] 
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D> D 

A COS [^(2 arcsin) + ^(2 arcsin) — 7r/4] 



rz> I 

A cos [^(2 arcsin) + / y^(-5;>)/(-A;>) 



7r/4] 



D. 



A COS [^(2 arcsin) - 2y^(-^>)V(-A;>)/3 



7r/4]. 



(81) 



"arcsin" denotes the arcsin function from eq. (|65p . Comparing 
eqs. dHU) and dZil), I find 



I have used the more accurate expression, eq. for the arcsin. 

The formula contains utt, rather than 2n7r, because c> can have 
either sign. Except for the sum and the unfamihar arcsin, the con- 
nection formula has the usual form. Again, the turning points are 
the two roots of — = which are closest to the center of the 
Mexican hat. (For more accuracy, use eq. and find the two 
roots of (7_ + — 2A = 0.) 

This is a good point to revisit the case briefly. If one 

thinks of the Mexican hat as an upside-down potential, then, be- 
cause n is large at A^ 0, one should expect the form eq. to 
work well, not badly. However,in the WKBJ approach, one always 
solves the equations several times: once away from turning points 
and once at each turning point. From eq. (f7T|) . 1/k blows up at 
a turning point which is also an SU(2) zero, since g_ has a square 
root zero there. The problem as A^ ^ is therefore at the turning 
points. They need a more careful discussion which I do not give 
here. 

If convenient for numerical purposes, I can replace the sum by 
an integral, as I did when computing S near turning points. The 
arcsin is monotonic and positive, so that (by the same arguments 
as those used to bound series by integrals in the Cauchy integral 
test for convergence) one can bound the sum above and below by 
two integrals with lower limits differing by one unit. These bounds 
should be fairly tight, if the integral extends over the entire classi- 
cally allowed region of order L. (I needed to do numerical work at 




2 arcsin J {g^ + g+ - 2A)/(4g_)dD = (n + l/2)7r. 




(82) 
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eq. (fTSj) . only because the integral was confined to a limited region 
and near turning points, where presumably the integral bounds are 
least restrictive). 

I attempted to obtain an analytic form for A, first replacing sum 
by integral in eq. (j82p . then carrying our the integral. When an 
integrand contains an arcsin of complicated argument, usually it is 
best to integrate by parts, which replaces the arcsin by a square root 
(also complicated, but not as bad as the arcsin). 

/ 2 arcsin - A)/2^_rfD = -A/ dD D d(g)^ /dD 

- [2ig-r^{g-y - m (83) 

(The integrate by parts surface term vanishes.) Eq. (jHHj) is reducible 
to elliptic integrals, but I believe this fact is of little use for extracting 
an analytic form for A. The arguments of the elliptic integrals will 
be functions of the zeros of g"^ — X"^, so that the A dependence will 
be implicit. 

As a numerical check on calculations of this section, I solved the 
exact difference equation on a computer for the case F = and 
Lx = Ly = L = 5. I obtained the eigenvalues 

A = 26.6, 19.6,13.3,7.99,3.59, (84) 

plus five more which were the negatives of these eigenvalues, plus one 
zero eigenvalue. I then used the FindRoot command on Mathemat- 
ica to solve the WKBJ quantization condition eq. (jH^ . FindRoot 
must be given n, plus two initial guesses as to the value of A; Find- 
Root then iterates, Newtonian style, until the program converges to 
a solution. I obtained 

A = 26.2,19.3,13.1,7.86. (85) 

The four values correspond to n = (largest) through n = 3 (small- 
est). I was unable to compute the n = 4 eigenvalue (the A^ — > 
eigenvalue). The FindRoot program would not converge. The agree- 
ment between eqs. (I84|) and (j85j) is surprisingly good for this small 
value of L. 
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